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We show that the decay of the inflaton field may be incomplete, while nevertheless successfully 
reheating the universe and leaving a stable remnant that accounts for the present dark matter abun¬ 
dance. We note, in particular, that since the mass of the inflaton decay products is field-dependent, 
one can construct models, endowed with an appropriate discrete symmetry, where inflaton decay is 
kinematically forbidden at late times and only occurs during the initial stages of held oscillations 
after inflation. We show that this is sufficient to ensure the transition to a radiation-dominated 
era and that inflaton particles typically thermalize in the process. They eventually decouple and 
freeze out, yielding a thermal dark matter relic. We discuss possible implementations of this generic 
mechanism within consistent cosmological and particle physics scenarios, for both single-held and 
hybrid inflation. 


I. INTRODUCTION 

Inflaton and dark matter candidates in particle physics 
models share several common features, both being typi¬ 
cally assumed to be weakly interacting and neutral fields. 
The inflaton scalar field must interact weakly with itself 
and other degrees of freedom in order to ensure the re¬ 
quired flatness of the associated scalar potential, which 
could be spoiled by large radiative corrections |l|. Sim¬ 
ilarly, dark matter particles should form a stable non- 
relativistic and non-luminous fluid at late times that ac¬ 
counts for the observed galaxy rotation curves [2] and 
the large scale structure in the universe as inferred from 
Cosmic Microwave Background a a and weak-lensing 
observations 0 fe] (for a recent review see 0)- Both 
inflation and dark matter are essential features in the 
modern cosmological paradigm and cannot be accounted 
for within the present framework of the Standard Model 
of particle physics. It is therefore interesting to consider 
the possibility that the same field accounts for both ac¬ 
celerated expansion in the early universe and the hidden 
matter component at late times. 

Scalar fields have the interesting property of mimicking 
fluids with different equations of state depending on the 
kinematical regime considered. For a homogeneous scalar 
field <f> with potential V (</>), we have: 

P<t>=\<i> 2 + V(<t>) , = (1.1) 

Hence, on the one hand a slowly-varying field, </> 2 /2 -C 
V(4>) acts as an effective cosmological constant, which is 
the regime typically considered in canonical inflationary 
models. On the other hand, a field oscillating about the 
minimum of its potential where U(</>) ~ m 2 </> 2 /2 behaves 
as non-relativistic matter, with (</> 2 /2) = —{V((/))) such 
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that p$ -C [§]. These two regimes will generically 
be present in inflationary potentials, which further sug¬ 
gests a common framework for inflation and dark matter. 
The main difficulty in achieving such a unified descrip¬ 
tion lies, however, in the fact that inflation must end with 
a transition to a radiation-dominated universe, in order 
to recover the standard “Big Bang” evolution at least 
before the freeze-out of light nuclear abundances takes 
place [9}. An efficient transfer of energy between the 
inflaton field and radiation generically requires the for¬ 
mer to decay into light degrees of freedom following the 
period of inflationary slow-roll even though other 
non-perturbative processes such as parametric resonance 
amplification could contribute significantly to the reheat¬ 
ing process [llj . A plausible alternative is also that non¬ 
equilibrium dissipative processes continuously source a 
radiation bath during inflation that smoothly takes over 
as the field exits the slow-roll regime, a possibility gener¬ 
ically known as warm inflation [12l I13| . 

Nevertheless, efficient reheating does not imply the 
complete decay of the inflaton field and the possibility 
that a stable remnant is left after inflation has been con¬ 
sidered in a few studies in the literature (l3 - 0 |. In this 
work, we propose a concrete realization of this generic 
idea in quantum field theory, where the decay of the infla¬ 
ton is truly incomplete, occurring only for a finite period 
after the end of the slow-roll regime. 

We describe a simple mechanism based on standard 
Yukawa coupling between the inflaton and fermion fields 
(and potentially their superpartners), endowed with an 
appropriate symmetry that protects the full decay of the 
inflaton field, and discuss the parametric regimes where 
inflaton decay into such fermions is incomplete. We also 
discuss the embedding of this generic mechanism in con¬ 
crete inflationary models and possible scenarios for the 
interactions between the fermion fields and the Standard 
Model degrees of freedom that allow for the presence of 
the latter in the post-inflationary thermal bath. We also 
show that the inflaton remnant is not necessarily in the 
form of a coherent condensate of bosonic particles and 
that, in particular, for not too small couplings the evap- 
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oration of this condensate is inevitable. In this case, the 
thermalized inflaton particles eventually decouple from 
the radiation bath and their abundance freezes out, yield¬ 
ing a thermal inflaton relic with properties similar to 
other WIMP candidates 0 - We suggestively denote 
this as the “WIMPlaton” scenario. 

In this work, we consider two dynamically distinct sce¬ 
narios. Firstly, we will consider a minimal model with a 
single dynamical scalar field that simultaneously drives 
inflation, reheats the universe through incomplete decay 
and leaves a stable non-relativistic remnant. Secondly, 
we consider a (supersymmetric) hybrid inflation model 
with an additional dynamical waterfall sector, which is 
responsible for ending inflation and reheating of the uni¬ 
verse. In this case, we show that despite reheating be¬ 
ing ensured by a different field, the inflaton field must 
decay in order for radiation to fully come to dominate 
the energy balance in the universe. As in the minimal 
model, this decay may nevertheless be incomplete and 
we discuss the parametric regimes in which the inflaton 
remnant constitutes a suitable dark matter candidate in 
the hybrid framework. 

This work is organized as follows. In the next section 
we discuss the minimal model of inflaton dark matter 
from incomplete decay, starting with its basic properties 
and dynamics and then discussing the possibility of con¬ 
densate evaporation, its embedding in a consistent infla¬ 
tionary model and different phenomenological possibili¬ 
ties for reheating the Standard Model degrees of freedom. 
In Section III we discuss these aspects in the hybrid re¬ 
alization of inflaton dark matter. Finally, in Section IV 
we summarize our main results and conclusions. 


II. MINIMAL MODEL 
A. Basic properties and dynamics 

The minimal model for inflaton dark matter considers 
a single dynamical (real) scalar field, the inflaton </>, with 
a potential energy V((f)) such that a period of slow-roll 
can occur for some field range. We take the inflaton held 
to be coupled to fermion fields ■)/’+ and -0 _ through stan¬ 

dard Yukawa terms and impose a discrete symmetiyQ [20| 
C 2 C Z 2 x S 2 on the Lagrangian such that the scalar in¬ 
flaton transforms under the Z 2 group as <fr —/ —</> and the 
fermions are simultaneously interchanged by the permu¬ 
tation symmetry ip + o if>-. This yields for the resulting 
Lagrangian density: 

£ = V{\(f)\) 

+ ^+(* 7 ^^ - mf)i> + + (* 7 ^ < 9/4 - m/)^_ 

— h(f)'if> + 'tl) + + h^-ifj- , (2.1) 


1 This subgroup contains only elements that transform simultane¬ 

ously under Z2 and S 2 . 


where we note that, as a result of the discrete symmetry, 
the two fermions have the same tree-level mass to/ but 
opposite Yukawa coupling to the inflaton held. The ac¬ 
tion of the discrete symmetry is restricted to the inhaton- 
fermion sector, such that all other helds, including the 
Standard Model helds, are invariant under this symme¬ 
try. 

This symmetry implies that the inflaton potential is 
an even function of the held and we assume that it is 
unbroken at the minimum of the potential, which must 
therefore lie at the origin, <f> = 0. As a result, since no 
other terms linear in cf> except for the above Yukawa terms 
are allowed by the discrete symmetry, the only possible 
decay channels of the inflaton at the minimum are (f> —> 
V’iV'i- If the inflaton mass at the minimum is given by 
to^ = V"(Q) > 0 and we require to/ > 771^/2, these 
decays will be kinematically forbidden for cf> = 0 and the 
inflaton will be stable when it settles at the minimum 
of its potential. Hence, the discrete symmetry and the 
mass hierarchy ensure the stability of the inflaton at late 
times, such that it will contribute to the present dark 
matter abundance. 

The two fermion fields are indistinguishable apart from 
their coupling to the inflaton, which leads to a mass split¬ 
ting away from the origin: 

m± = |to/ ± h(f>\ . ( 2 - 2 ) 

This implies that, although inflaton decay into these 
fermions is forbidden at late times as the field approaches 
the origin, that need not be the case if the amplitude of 
field oscillations after inflation is sufficiently large. In 
particular, decay will be allowed for field values satisfy¬ 
ing: 


|to/ ± h(f >| < to^/2 . (2-3) 


This implies, in particular, that decay is kinematically al¬ 
lowed for field amplitudes \<f>\ > rrif/h. The partial decay 
widths associated to the two fermionic decay channels are 
then given by: 


F h 2 

r± = 



(2.4) 


with r 0 = F+ + F_. Note that due to the opposite sign of 
the Yukawa couplings, the inflaton will alternately decay 
into each fermion species as it oscillates between negative 
and positive values. The inflaton equation of motion as 
it oscillates about the minimum of its potential is then 
given by: 

(f> + 2>H(f) + Y^cf) + rn^cf) = 0 , (2.5) 

which, upon multiplying by (f> is equivalent to: 

p<j> + 37L(p0 + P<t>) = -F </<(> 1 2 . (2.6) 

The term on the right hand side gives the rate at which 
energy density is transferred from the inflaton field into 
the fermions (f’±. Let us assume that the fermions quickly 
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thermalize, an assumption that wc will check a posteriori , 
in the process exciting g * relativistic degrees of freedom 
and forming a radiation bath at temperature T, with 
energy density pr = (7t 2 /30 )g*T 4 . Note that since the 
fermion masses oscillate due to the varying inflaton field, 
they only contribute periodically to the number of rela¬ 
tivistic degrees of freedom. These may also include the 
inflaton for T > and other species such as the Stan¬ 
dard Model particles. Note that the latter must be ex¬ 
cited before the cosmological synthesis of light nuclear 
elements takes place, as we discuss in more detail in sec¬ 
tion IID. For simplicity, we consider a fixed value of < 7 *, 
which is not a bad approximation since, as we will show 
later on, our results exhibit only a mild dependence on 
this parameter. 

Energy conservation then implies that the energy lost 
by the inflaton field in Eq. (12.61) is gained by the radiation 
bath, which then follows the dynamical equation: 

PR + AH PR = T^ 2 . (2.7) 


Since both the inflaton and the radiation contribute to 
the energy density in the universe, we may write the 
Friedmann equation as: 


_ pp + pr _ \'<\> 2 + h™ 2 ^ 2 + PR 
3 Mf, ~ 1™ 2 


Eqs. (12.51) . (12.71) and (12.81) then form a complete set of 
differential equations that can be solved for given choices 
of the parameters (m^, mf,h) and initial conditions. The 
non-linearity of the equations precludes a complete an¬ 
alytical description of the problem, although as we dis¬ 
cuss below a few simple analytical estimates can be per¬ 
formed, but it is straightforward to integrate them nu¬ 
merically. It is useful to write the fermion tree-level mass 
as: 

m f = ^(l + S), (2.9) 


such that the decay is kinematically forbidden (allowed) 
at late times for <5 > 0 (< 0 ). 

Figure [T| shows an example of the results obtained by 
numerically solving the inflaton-radiation equations for 
m<f) = 10 _ 3 Mp, S = 0.1 and h = 1. In this example we 
take as initial conditions at t = ti = 10 _2 to7 , (f>(ti ) = 

Mp, <p(ti) = 0 and PR{ti) = 10 ~ 16 p ( j,(t i ), although the 
results are essentially unchanged as long as the field value 
is sufficiently large before oscillations begin for the decay 
to be kinematically allowed as discussed above and the 
radiation energy density is initially negligible, as should 
be expected after 50-60 e-folds of inflation. 

The plot in Figurc[T](a) illustrates the main dynamical 
features that are generically obtained. The inflaton field 
begins to oscillate about the origin with frequency 
for t ~ m ^) 1 and behaves initially as cold dark matter, 
p^ oc a -3 , being the dominant energy component such 
that a oc t 2 ^ 3 and H = 2/3 1. While decay is blocked be¬ 
fore the onset of oscillations, since m± ~ h\<j)\ m^, it 


(a) 



(b) 



(c) 
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FIG. 1. Results of the numerical integration of the inflaton- 
radiation dynamical equations for = 10~ 3 Mp, 5 = 0.1 
and h = 1, showing the time evolution of (a) the inflaton 
(solid blue curve) and radiation (solid red curve) energy den¬ 
sities; (b) the inflaton-to-entropy ratio; and (c) the inflaton 
decay width (solid blue curve) compared to the Hubble pa¬ 
rameter (solid green curve). The blue dashed curves in (a) 
and (c) give the evolution of the inflaton energy density in 
the absence of decay and the maximum value of the decay 
width, respectively. All quantities are given in Planck units 
such that Mp = 1. 
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becomes kinematically allowed as soon as the field goes 
through the origin. Since in this example the ampli¬ 
tude of the field oscillations largely exceeds the tree-level 
masses, with 3> m/, decay occurs for two narrow 
field ranges close to (and on both sides of) the origin. 
The decay width then corresponds to the series of peri¬ 
odic narrow peaks shown in Figure [Q (c), with maximum 
value T™ ax = (h 2 / 87 r)m</. While initially T™ ax <C 3 H 
as illustrated in Figure Q] (c), such that the inflaton’s 
energy density remains essentially unaffected by the de¬ 
cay into fermions, the source term in the radiation equa¬ 
tion quickly becomes significant, leading to a jump in the 
value of the radiation energy density. The latter remains 
approximately constant until the inflaton’s energy den¬ 
sity is sufficiently redshifted. When they become com¬ 
parable in magnitude, the inflaton effectively decays and 
radiation takes over as the dominant component. As the 
field amplitude decreases, the maximum decay width be¬ 
comes progressively smaller until decay is finally blocked. 
The oscillating inflaton then becomes stable and once 
more behaves as cold dark matter, eventually taking over 
the radiation as the dominant component at later times. 

A peculiar dynamical feature of the evolution is the 
approximate constancy of the radiation energy density 
achieved just after the first few oscillations. This is in¬ 
herent to the fact that inflaton decay occurs in short 
bursts in each oscillation, which does not occur if the de¬ 
cay were allowed for all field values. Before the effect of 
decay into fermions becomes significant, for t > m ^ 1 the 
inflaton behaves as a damped harmonic oscillator with 

(pit) ~ $(f) sin(m 0 t + a<f,) , $(f) = 

( 2 . 10 ) 

where a^ is a phase depending on the initial field and 
velocity values. For /i<F to/, one can easily see that 
decay into each fermion is allowed during a short period 
Td ~ (/id >) -1 < 27t/to/, as the held goes through the ori¬ 
gin, which occurs twice every oscillation period. Since 
the average decay width in this interval can be taken as 
T ™ ax /2 and the held velocity is (p ~ >, every half pe¬ 

riod the radiation energy density increases due to inflaton 
decay by an amount: 

pmarr 7 

a p d R Cav ^ -y- (TO^$) 2 (2r d ) ~ — ml$ , (2.11) 

where we have taken into account the decays into both 
ip+ and ip-. When decay is forbidden, radiation simply 
redshifts with expansion, which counteracts the enhance¬ 
ment due to decay by an amount: 

A p* ubble ~ -AHp R (n/m^) . (2.12) 

Since H = 2/3f = (to<//-\/6Mp)$, it is easy to see that 
the amount of radiation produced by inflaton decay can 
be balanced exactly by Hubble expansion to yield a con¬ 
stant energy density. Equating A p d ^ cay = — Ap Pubble 


l~8 M P 1 
V 3 nip ~t 


then gives: 

PR - ^2 m l Mp ’ ( 2 - 13 ) 

which is in very good agreement with the numerical sim¬ 
ulations. Noting that, from Eq. (mi), p R ~ 0 implies 
4 Hp R ~ Tp<p 2 ~ r ppp, we see that that ~ 3 H for 
p R ~ Ptf) , so that as observed numerically the inflaton en¬ 
ergy density is only reduced signihcantly when it becomes 
comparable to the radiation energy density. 

From Eq. (12.131) we can easily determine the associated 
temperature, which remains approximately constant up 
to inflaton-radiation equality and thus constitutes the re¬ 
heating temperature: 


~ 2.7 x 10V 17 V /4 3/4 GeV. (2.14) 

From this we conclude, in particular, that the reheating 
temperature is generically larger than the inflaton mass 
for: 

mp < 5.6 x 10 16 — GeV . (2.15) 

9* 


After inflaton-radiation equality, the field decays ex¬ 
ponentially fast until its amplitude drops sufficiently for 
decay into fermions to become inefficient. Numerically, 
we observe that this occurs before the decay becomes 
kinematically forbidden, for $ < rrif/h , corresponding 
roughly to when the (maximum) decay width becomes 
less than the Hubble rate. Since afterwards the field sta¬ 
bilizes and behaves as non-relativistic matter, it is useful 
to compare its number density np = pp/rrip oc a~ 3 with 
the radiation entropy density: 


2 tt 2 

S = ~45 9 * T3 ’ (2 ‘ 16) 

which redshifts in a similar way up to changes in the 
number of relativistic degrees of freedom, with T oc a -1 . 
As illustrated in Figured (b), the ratio np/s stabilizes af¬ 
ter decay becomes inefficient, and numerically we obtain 
the following expression: 

^ ~ 0.5 ff “ 1/4 h- 3 - 6 (mp/Mp) 1 ' 02 f{6) , (2.17) 
where for 0 < 5 < 1: 

f(5) = 1 + 4.8v^ + 0.5(5 . (2.18) 


If no other processes change the inflaton particle number 
density in the oscillating field, as we investigate below 
in more detail, this ratio will remain constant until the 
present day. Assuming the oscillating inflaton field ac¬ 
counts for all the dark matter in the universe, the present 
inflaton-to-entropy ratio is given by: 

n <f >o _ 3 Tp H c q g*o ^ 2Q—28 Mp p ig) 

Sq 4 mp Q R q <?* s o mp 
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where To = 2.73 K is the present CMB temperature, 
while Yl c0 h 2 = 0.12 and Ylp 0 h 2 ~ 4.17 x 10 -5 are the 
present abundances of cold dark matter and radiation, 
respectively Q. We have also taken into account the 
present difference between the number of relativistic de¬ 
grees of freedom contributing to the radiation and en¬ 
tropy densities, g *o — 3.36, g * s o ~ 3.91. Equating (12.171) 
and (12.191) . we obtain for the inflaton mass: 

rritf, ~ 72 g° 12 h 1/78 f(S) 1A9 TeV . (2.20) 

For fermion masses above the kinematical limit and of 
the order of the inflaton mass, 0 < <5 < 1, and taking 

= 10 — 100, we conclude from Eq. (12.141) that the re¬ 
heating temperature is above 100 MeV for Yukawa cou¬ 
plings h > ICY 6 — 10 -5 , corresponding to inflaton masses 
m</ > 10 keV. For larger (and arguably more natural) 
couplings h > 10~ 3 — 10~ 2 , the inflaton may account 
for the dark matter in the universe for masses in the 
GeV—TeV range, similarly to the mass range obtained 
for WIMP-like thermal relics. Qualitatively, it is easy to 
understand the parametric dependence of the required 
inflaton mass. The (incomplete) decay into fermions is 
more efficient for larger couplings, which affect both the 
overall value of the decay width and the effective fermion 
masses, and lighter fermions, yielding a smaller inflaton- 
to-entropy ratio at late times and thus allowing for larger 
inflaton masses to match the present dark matter abun¬ 
dance. On the other hand, smaller couplings and heavier 
fermions lead to a larger inflaton abundance, which may 
overdose the universe unless the inflaton is sufficiently 
light. 

We note also that there is no gain in considering finely- 
tuned fermion masses, i.e. S « 0, since even though 
decay is kinematically allowed for longer as d decreases, 
decay becomes inefficient when T/, ~ H, a condition that 
becomes independent of the tree-level fermion mass in 
this limit. Although we have restricted the numerical 
analysis to values of 6 < 1, we expect incomplete de¬ 
cay to be efficient as long as T^ > H for inflaton field 
values satisfying the kinematical condition in Eq. m - 
In particular, since the decay width takes its maximum 
value Y™ ax = (h 2 /8n)m t j > for h\<f>\ = m/, i.e. when the 
fermions are effectively massless, and H ~ m^ipl/Mp 
during inflaton-domination, we conclude that decay will 
be efficient for: 

to/ < h 3 M P , (2.21) 

which allows for quite large fermion masses if the Yukawa 
coupling is not too suppressed. 

B. Condensate evaporation: the WIMPlaton 
scenario 

Although the imposed discrete symmetry protects the 
inflaton from decaying into any other particles except for 
the fermions the above analysis neglects the effects 
of additional interactions induced by the Yukawa terms 


in Eq. EH) and which may play an important role as we 
discuss below. 

The classical inflaton field corresponds to a collective 
state of zero-momentum scalar bosons, assuming that 
no large field inhomogeneities are formed at the end of 
the slow-roll inflationary regime. Inflaton particles in 
this condensate can interact with the fermions that re¬ 
sult from its decay and, in particular, these fermions can 
scatter some of the bosons out of the condensate and 
promote them to higher-momentum states that become 
part of the thermal bath. These correspond to scatter¬ 
ing processes ip±{4>) —> ip±4>, where we denote by (</>) 
and 4> scalar particles in the zero-momentum condensate 
and in higher-momentum modes, respectively, and which 
are mediated through both s- and f-channel fermion ex¬ 
change. Moreover, these processes may occur as soon as 
the field begins oscillating and decay into ip± becomes 
kinematically allowed, potentially leading to the evapo¬ 
ration of the condensate and the transfer of the inflaton 
particle number into the thermal bath. 

As we have seen earlier, soon after the onset of infla¬ 
ton oscillations, the temperature of the thermal bath rises 
sharply to a value that remains approximately constant 
until inflaton-radiation equality and that corresponds to 
the reheating temperature in Eq. ( 12 . 141 ) . In particular, 
this temperature exceeds the inflaton and fermion (tree- 
level) masses in the parameter space region that yields 
the present-day dark matter abundance. Further assum¬ 
ing that local thermal equilibrium is quickly achieved, as 
we check below, we may then take the phase-space dis¬ 
tributions for inflaton and fermion species in the thermal 
bath to be the relativistic Bose-Einstein and Fermi-Dirac 
distributions, respectively. Taking into account the above 
scatterings and the inverse processes, the net condensate 
evaporation rate is given by [2l| : 

t -°- =£ / n (Sjfe 1 p* +P2-P3-P,) 

X|M| 2 lhk( 1 + /„)( 1 - h) - hi 4(1 + /i)(l - /„)] , 

( 2 . 22 ) 

where M is the scattering amplitude for {(/>) (pi)ip±(P 2 ) ^ 
<p(p 3 )ip±{ pa) and /) the corresponding phase-space distri¬ 
bution factors. Since the condensate is inherently char¬ 
acterized by large occupation numbers /i > 1, we obtain 
to leading order for T to,/, to±: 

TeVap ~ 1^3 ( 1 + l0S (^)) T ’ < 2 ' 23 > 

where we have taken into account the contribution of 
both fermion species. Note that this is only valid when 
the fermions are relativistic, while for non-relativistic 
fermions in local thermal equilibrium the evaporation 
rate is exponentially suppressed. This means that during 
the first few oscillations before the inflaton amplitude is 
significantly reduced, condensate evaporation only occurs 
during the short periods where decay is also kinematically 
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allowed. As we have concluded above, for the parameter 
values yielding the present dark matter abundance we 
obtain Tp> > and since m± ~ to/ after reheat¬ 

ing the above expression holds until either the inflaton or 
the fermions become non-relativistic. Since in the radi¬ 
ation era H ~ (ir/y/90)gl 2 T 2 /Mp, condensate evapora¬ 
tion becomes progressively more efficient as the temper¬ 
ature drops. We may then determine a lower bound on 
the Yukawa coupling such that condensate evaporation 
is inefficient for T > to/,, to/. For comparable masses, we 
obtain: 


r 

evap 

H 


*L-(£±Y 1/2 

47T 3 \ 10/ \ TO/, ) 


< 1 . (2.24) 


Using Eq, (12.201) for the inflaton mass, we then obtain for 
<5 = 1: 


h < 10 ~ 5 g° 28 . (2.25) 

As we had seen above, this is in tension with the lower 
bound on the coupling required for a reheating tempera¬ 
ture above 100 MeV, such that initial conditions for Big 
Bang Nucleosynthesis (BBN) are already in place after 
inflaton decay JUJ. This conclusion is essentially com¬ 
mon to all fermion mass values to/ > to/,/ 2. Although 
a more detailed analysis of the Boltzmann equation de¬ 
termining the evolution of the inflaton condensate may 
be required, this estimate indicates that in the physically 
interesting parameter range, where the condensate could 
account for the present dark matter abundance while sat¬ 
isfying the BBN constraint, condensate evaporation is 
most likely inevitable. 

Since evaporation simply transfers the zero-momentum 
condensate particles into excited states, the conclusion 
above does not imply that inflaton particles cannot ac¬ 
count for dark matter. One can easily check that all 4- 
body processes involving relativistic inflaton and fermion 
particles, including scatterings and annihilations, occur 
at a rate comparable to the evaporation rate obtained in 
Eq. (12.231) . Hence, evaporation of the condensate should 
lead to a bath of fermions and inflaton particles (as well 
as other species) in local thermal equilibrium. In fact, 
this may even occur just after the onset of inflaton os¬ 
cillations, for H < to/,, when the temperature reaches 
an approximately constant value T > m^^mf as seen 
above, for h < 0(1) couplings. The condensate’s energy 
is in this case quickly transferred to the thermal bath, 
increasing its temperature to a maximum value: 

Tr 1x = ( 5 ) ' 

- 8-5 x 10 10 ff r 1/4 GeV . (2.26) 

Note that this was obtained for H = to/,, which con¬ 
stitutes only an upper bound since fermion particles are 
only produced for H < to./, as shown by our numerical 
simulations. The true reheating temperature will then lie 


between the value obtained in Eq. (12.141) and this max¬ 
imum value if condensate evaporation occurs before it 
decays significantly. The reheating temperature is given 
by Eq. (12.141) if condensate evaporation only occurs in the 
radiation era, which corresponds to Yukawa couplings 
h < 10” 4 for <5 < 1 and g* = 10 — 100. Note that in 
this parametric regime thermalization of the fermions in 
the plasma can only be efficient if their interactions with 
other degrees of freedom are stronger than those induced 
by the Yukawa terms, while for h > 1CU 4 the latter occur 
sufficiently fast to maintain local thermal equilibrium, as 
we analyze in more detail in section IID. 

After reheating, inflaton and fermion particles will be 
kept in local thermal equilibrium by annihilation and 
elastic scattering processes. Once these become inef¬ 
ficient, the abundance of inflaton particles will freeze 
out, as for other conventional WIMP dark matter candi¬ 
dates. Assuming this occurs when both the inflaton and 
the fermions are non-relativistic, the relevant (fermion 
t-channel) annihilation cross section is given by: 


& cf>(f> — 


h A 

8nm 2 ’ 


(2.27) 


which is independent of the fermion mass in this limit. 
Following the standard calculation for the thermal relic 
abundance of a decoupled non-relativistic species, we ob¬ 
tain for the inflaton mass: 



(2.28) 

where g„p denotes the number of relativistic degrees of 
freedom at freeze-out and xf = ru^/Tp, with Tp denot¬ 
ing the freeze-out temperature. This is somewhat smaller 
than the mass values obtained assuming the oscillating 
inflaton condensate survives until the present day, al¬ 
though in a comparable range and exhibiting a similar 
dependence on the Yukawa coupling. 

This then gives us a more realistic dynamical picture 
of what we suggestively denote as the “WIMPlaton” sce¬ 
nario. After inflation, the scalar inflaton begins oscillat¬ 
ing about the minimum of its potential, decaying into 
fermions in short bursts every oscillation. These may 
thermalize and excite other degrees of freedom in the 
plasma, and scatter off the inflaton particles in the con¬ 
densate, leading to its evaporation. Both decay and evap¬ 
oration increase the relative abundance of radiation and 
decrease the amplitude of oscillations, until eventually ra¬ 
diation becomes dominant and inflaton decay is no longer 
kinematically allowed. The stable inflaton particles re¬ 
main in thermal equilibrium until the temperature drops 
below their mass and they decouple from the plasma, 
their frozen abundance yielding the inferred dark matter 
component of our present universe. 
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C. Embedding in a consistent inflationary model 


As we have concluded from the analysis above, the in- 
flaton field can account for dark matter in the universe 
at late times for masses below or around the TeV scale in 
the WIMPlaton scenario, with a similar mass range ob¬ 
tained assuming the inflaton condensate does not evap¬ 
orate. This implies that the inflaton potential cannot 
be given solely by a quadratic mass term, since the am¬ 
plitude of CMB temperature anisotropies would yield in 
this case: 


fn<l> — 


V 6it2A x 


N e 


Mp ~ 1.4 x 10 


13 



GeV , 

(2.29) 


which for 50 — 60 e-folds of slow-roll inflation largely ex¬ 
ceeds the TeV scale. However, our analysis assumed only 
that is the inflaton mass as it oscillates about the 
minimum of the potential at the origin, while the effec¬ 
tive inflaton mass can be much larger if slow-roll occurs 
for significantly larger field values where self-interactions 
become important. For example, the discrete Z 2 symme¬ 
try allows for quartic self-interactions such that: 

V{<t>) = > (2-30) 

with the quartic term dominating for \cf>\ > 

Since it increases the inflaton mass, the quartic self¬ 
interaction would make the decay into fermions more effi¬ 
cient and speed up the reheating process. Our analysis is 
nevertheless valid if the quartic term is sub-dominant for 
the field values \(j)\ ~ 'mf/h at which decay into fermions 
occurs, which requires: 


h> 


,/A (ati) = 3 x 10 - s (mi) , 

V 12 } \m 


(2.31) 


where we have used A ~ 10 -14 as imposed by the COBE 
normalization for inflation with a quartic potential. This 
is easily satisfied if the fermions are not much heavier 
than the inflaton given the more stringent bounds on 
the Yukawa coupling discussed earlier. The addition of 
a quartic term is, however, not sufficient to produce a 
consistent model of inflation, since it predicts a too red- 
tilted spectrum for curvature perturbations and a tensor- 
to-scalar ratio already outside the bounds obtained by 
Planck |23| and BICEP2 [U|. 

A consistent spectrum may be achieved in warm infla¬ 
tion scenarios [l2| , where the coupling between the infla¬ 
ton and other sectors plays a significant role during in¬ 
flation. For example, renormalizable couplings between 
the inflaton and additional fermionic and bosonic species 
may lead to dissipation of the inflaton’s energy into ra¬ 
diation during the slow-roll period, leading to an obser- 
vationally consistent spectrum if inflaton particle states 
thermalize with the other relativistic degrees of freedom 
[25l | . Since this naturally leads to radiation becoming the 
dominant component at the end of the slow-roll regime, 


the post-inflationary evolution will necessarily differ from 
the dynamical picture discussed in this work. It would 
nevertheless be interesting to investigate whether the in¬ 
flaton can act as a stable dark matter candidate at late 
times in this context, although this lies outside the scope 
of the present discussion. 

Another interesting possibility is the inclusion of a non- 
minimal coupling to the gravitational sector, in particular 
a coupling between the inflaton and Ricci scalar of the 
form which is compatible with the discrete Z 2 sym¬ 

metry. The resulting inflationary scenario yields a per¬ 
turbation spectrum that smoothly interpolates between 
the minimal quartic model and the Starobinsky model 
as the non-minimal coupling constant increases. On the 
one hand, the latter is characterized by a low tensor- 
to-scalar ratio and a spectral index n a = 0.96 — 0.97 in 
agreement with the Planck results; on the other hand, 
a small non-minimal coupling constant is preferred to 
obtain a non-negligible tensor-to-scalar ratio . We refer 
the reader to [26|, (27| for a more detailed discussion of 
these scenarios, since here we are mainly interested in 
the post-inflationary dynamics. The effect of the non- 
minimal coupling on the effective scalar potential in the 
Einstein frame becomes negligible for £</> 2 /Mp < 1, such 
that consistency of our analysis implies 


h> \/e 


M P ’ 


(2.32) 


which is generically less stringent than Eq. (12.311) for 
masses in the TeV range and £ < 10 16 . 

Finally, an important aspect in embedding the in¬ 
teractions in Eq. m within a consistent inflationary 
model is the fact that the discrete Z 2 x S 2 symmetry 
does not protect the scalar potential from radiative cor¬ 
rections. In particular, the Yukawa interactions induce 
loop-corrections of the Coleman-Weinberg form, which 
for large inflaton field values take the leading form: 

/i 4 </> 4 (i ^(h 2 (j) 2 \ 3' 


*V f «- 


16-7T 2 


log 


(2.33) 


and therefore induce an effective quartic term in the 
potential. The effect of this term does not necessarily 
spoil the predictions of the non-minimally coupled quar¬ 
tic model, as discussed in [28|, although one must en¬ 
sure that the observed normalization of the perturbation 
spectrum is obtained. While for ( < 1 the effective quar¬ 
tic coupling must have approximately the same value as 
in the minimally coupled case, which requires h < 10 , 

significantly larger values can be accommodated for large 
non-minimal couplings. 

Radiative corrections can, however, be significantly re¬ 
duced in supersymmetric scenarios, and a supersymmet¬ 
ric version of the model in Eq. CU with a C 2 C Z 2 x 5 2 
symmetry can be easily obtained by considering a super¬ 
potential of the form: 


W = (Y 2 — Y 2 ) + (Y 2 + Y 2 ) + 


rritj, 




-<h 4 Z 


(2.34) 









where the inflaton and fermions ij)± are embedded within 
the chiral superfields $ and Y ±, respectively, and the aux¬ 
iliary superfield Z induces the quartic term in the infla¬ 
ton potential (noting that the discrete symmetry forbids 
cubic inflaton terms in the superpotential). Supersym¬ 
metry then cancels the leading contributions of scalars 
and bosons to the 1-loop Coleman-Weinberg potential, 
which becomes: 


A Vsusy ~ • (2.35) 

This contribution is thus necessarily smaller than the 
tree-level potential V((j>) ~ A 2 10| 4 /4 for h < 1, therefore 
avoiding the generation of large effective self-interactions 
during inflation. 

Besides the Yukawa terms considered so far, the super- 
symmetric model also yields scalar interactions between 
the inflaton and the scalar partners y± of the fermions 
ip ±, which apart from SUSY splittings that vanish at 
the origin have the same mass m± = \mf + h<f >|. Trilin- 
ear terms in the scalar potential also lead to the decay 
(j> —> y±y ±, with analogous kinematics and comparable 
widths to the fermionic decay channels, therefore yield¬ 
ing a similar incomplete decay of the inflaton as analyzed 
above. We note that the incomplete decay dynamics can 
be fully described in terms of scalar fields and is there¬ 
fore not exclusive of fermion Yukawa couplings, although 
the required form of the scalar mass terms is more nat¬ 
urally motivated within a supersymmetric context. We 
also note the existence of quartic terms in the scalar po¬ 
tential which induce the 3-body decay 4> —> zy±y ±, where 
z is the scalar component of the Z chiral multiplet. The 
associated couplings have opposite signs for y± and pro¬ 
portional to A, so that they are typically sub-dominant 
with respect to the 2-body decays and they are also kine¬ 
matically forbidden at late times, so that they do not 
affect our earlier conclusions. 


D. Reheating the Standard Model 

We have so far assumed that the fermions (or scalars as 
discussed above) resulting from inflaton decay thermalizc 
and excite other degrees of freedom, and in particular it 
is crucial that Standard Model particles are generated in 
the thermal bath at temperatures above ~100 MeV so 
that BBN may occur following the standard freeze-out 
dynamics of light nuclear abundances. 

As we have briefly discussed above, the fermions them¬ 
selves cannot be treated as fully relativistic degrees of 
freedom before the inflaton decays sufficiently, since their 
mass varies between small and large values as the infla¬ 
ton oscillates about the origin. For the short periods 
when they are light and decay is allowed, thermaliza- 
tion through Yukawa interactions can be quite efficient. 
For example, fermion-fermion scatterings through infla¬ 


ton s—channel exchange occurs at a rate: 


r 


ifjijj — 


9C(3) 

16-7T 3 


h A T . 


(2.36) 


Since after the onset of oscillations H < and tak¬ 
ing the temperature value obtained in Eq. (12.141) from 
inflaton decay, we have: 



> 


F 

me/, 


~ 60 S 7 1/ V ,/4 (iTfgy) 


- 1/4 


,(2.37) 


such that for 0(1) couplings Yukawa interactions can lead 
to thermalization from the start of oscillations. Recall 
that, as seen above, evaporation of the condensate oc¬ 
curs at a comparable rate, so that as the temperature 
rises due to evaporation both processes become progres¬ 
sively more efficient. If the Yukawa couplings are more 
suppressed, one may also envisage scenarios where for 
example the fermions are charged under a gauge sym¬ 
metry that is unbroken at the relevant temperatures and 
thermalization occurs through gauge boson exchange for 
sufficiently strong couplings. 

Outside the inflaton field range for which the fermions 
are effectively light, their number density will become 
Boltzmann-suppressed, and any interactions will neces¬ 
sarily become inefficient in keeping the fermions in local 
thermal equilibrium. We thus expect them to decouple 
for most of the oscillation period, transferring their en¬ 
tropy into light degrees of freedom such as the inflaton it¬ 
self or e.g. gauge bosons. If interactions occur faster than 
the inflaton oscillation rate ro^, as given in Eq. (12.371) for 
the Yukawa scattering processes, fermions will drop in 
and out of local thermal equilibrium as they oscillate be¬ 
tween the relativistic and non-relativistic regimes. This 
will lead to an oscillating g *, but as mentioned at the 
start of our discussion this is of little consequence since 
our results exhibit only a mild dependence on this pa¬ 
rameter. 

Since the inflaton is typically taken as a gauge singlet, 
the structure of the Yukawa interactions implies that if>± 
are non-chiral fermions (either Dirac or Majorana), as op¬ 
posed to the known SM fermions, being thus unlikely that 
they are explicitly charged under the SM gauge group. 
There are nevertheless several possibilities for exciting 
the SM degrees of freedom in the plasma either before or 
after radiation comes to dominate. 

The simplest case is perhaps that of unstable i/j± 
fermions decaying into a light scalar and a light fermion, 
for which the decay width is: 

h 2 f 

r v>± = ’ ( 2 - 3 8) 

where hf denotes the associated coupling constant. Note 
that this is computed in the fermions’ rest frame, whereas 
in the plasma’s frame an additional time dilation factor 
m±/T suppresses the decay for relativistic fermions. For 
a significant part of the energy density in radiation to be 
transferred into these degrees of freedom, the fermions 
must decay before they become non-relativistic. Thus, 
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requiring T^, ± > H for T > to/ in the radiation era 
(where m± ~ to/), we obtain the following bound on the 
coupling: 

6/>8xlO- s 4 4 ( I ^y 1/2 , (2.39) 

where g*/ is the number of relativistic degrees of free¬ 
dom at T = to/. In addition, for this to happen be¬ 
fore BBN we require to/ > 100 MeV. If the inflaton 
and fermion masses are comparable, this corresponds to 
h > 0.01 according to Eq. (12.281) . These light degrees 
of freedom may correspond to SM particles if, for ex¬ 
ample, the fermions coupled to the inflaton field corre¬ 
spond to a pair of degenerate sterile neutrinos, which are 
singlets under the SM gauge group and may decay into 
a Higgs-lepton pair through Yukawa terms of the form 
hjHlip±. Note that this requires to/ > tuh = 125 GeV 
[29| and hence < 250 GeV, which is compatible with 
the present dark matter abundance for couplings h < 0.7 
from Eq. (12.281) . 

If the fermions are stable, another possibility for re¬ 
heating the SM degrees of freedom is through efficient 
annihilation. A possible scenario is for the fermions to be 
charged under a hidden U(l)x gauge group, which may 
mediate fermion scatterings and thus improve the ther- 
malization efficiency. This U(l)x hidden photon may 
be kinetically mixed with the SM photon or hypercharge 
gauge boson through a term of the form F^F^ V , 
which may be generated radiatively if there are fields 
charged under both gauge groups or simply via gravita¬ 
tional interactions, as happens e.g. in string theory. Di- 
agonalization of the gauge kinetic terms then induces a 
small electric charge for the fermions i/j± , such that they 
may annihilate into SM charged particles via s-channel 
photon exchange, ip±ijj± — > 7 — > qq , ll ■ At high tempera¬ 
ture the annihilation cross section is given by the Thom¬ 
son scattering formula and the corresponding interaction 
rate for relativistic species is then given by: 

T th * ^-e 2 a 2 N ch T , (2.40) 

7r 


where e is the “mili-charge” of the fermions ?/±, a is the 
fine-structure constant and N c h is the effective number 
of charged species in the final state, which for the full 
SM is N c h = 20/3. In the radiation era, since FI oc T 2 
annihilation becomes more efficient at smaller tempera¬ 
tures. Once T < to/, however, the fermion abundance 
becomes Boltzmann-suppressed and annihilations can no 
longer be efficient. Thus, requiring that SM species are 
excited before ip± become non-relativistic, we obtain the 
following bound on the mili-charge: 


e > 5 x 1CT 7 





(2.41) 


from direct collider searches (including the LHC) yielding 
e < 0.1 for 1 GeV < to/ < few x 100 GeV and indirect 
bounds from the CMB anisotropy spectrum, which yield 
e < 10 -4 for fewx 100 GeV < rrif < fewxTeV based on 
the effect of mili-charged particles on the baryon-photon 
oscillations (for a gauge coupling gx = 0.1). Note that 
for masses above the TeV range, mili-charged particles 
may give a too large contribution to the dark matter 
abundance in the universe, and more stringent bounds 
on e apply in this case (see jloj and references therein). 
This thus constitutes a promising scenario with potential 
for experimental probing in the near future. 

The discrete C 2 C Z 2 x £ 2 symmetry protects the in¬ 
flaton from decaying at late times, thus constituting a 
viable dark matter candidate. One can consider, how¬ 
ever, scenarios where this symmetry is broken and the 
inflaton is only meta-stable, with a lifetime larger than 
the age of the universe, to ~ 14 Gyrs. Note that inter¬ 
actions between the fermions and other light degrees 
of freedom as in the scenarios outlined above can induce 
the decay of the inflaton through radiative effects or pro¬ 
cesses mediated by off-shell fermions. A few examples 
of these processes are illustrated in Figure [2] and in all 
cases the contribution of ip + and , 0 _ cancels if the dis¬ 
crete symmetry is exact. 

(a) (b) (c) 




FIG. 2. Feynman diagrams for the 2-body decay of the in¬ 
flaton into (a) gauge bosons, (b) light scalars and (c) light 
fermions, induced at the 1-loop level by gauge and Yukawa 
interactions of the tp± fermions. In (d) we also show the 4- 
body decay of the inflaton induced by the exchange of virtual 
ip± modes with Yukawa interactions with other light species. 
For clarity, all light fields are represented by green lines. 

As an example, we have considered the case where if>± 
are unstable, decaying into a light fermion and scalar, 
which induce diagrams (b)-(d) in Figure [2] Since these 
processes have comparable magnitudes, we have com¬ 
puted the inflaton decay width for the 1 -loop process 
(b), where it decays into two light scalars. For concrete¬ 
ness we consider the case where ip± are slightly non¬ 
degenerate with m_ = to/ and m+ = to/( 1 + A), for 
A <C 1. This gives for to/, = 2to/ (5 = 0): 


This bound is not very stringent for fermion masses in 
the GeV-TeV range, where the main constraints come 


r (6) = h2h ) 

^ 64-7T 5 


to</A 2 


(2.42) 
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which is more suppressed for larger values of to/. 

For 171^ yielding the correct relic abundance to account 
for dark matter and using the lower bound on hf obtained 
in Eq. (12.391) . we obtain the following upper bound on the 
inflaton lifetime: 


< 15.7 g*f 


f ™<t> y 4 1 

< A ^ 

VlOGeV/ 1 

v 0.01 / 


Gyrs ,(2-43) 


where we have taken the reference values for the present 
dark matter abundance and freeze-out parameters in 
Eq. (12.2811 . Hence, we conclude that significant violations 
of the discrete symmetry can still yield a sufficiently long- 
lived inflaton for 100 MeV < to/, < 10 GeV, such that 
the fermions i/j± decay before BBN and while they are 
still relativistic. For heavier inflaton particles, values of 
A below the percent level are required, signaling that the 
discrete symmetry most hold to a high degree of accuracy 
in this regime. 

It is thus clear from the examples above that if the in¬ 
flaton field can only decay to the fermions for a finite 
period following the end of the slow-roll regime, becom¬ 
ing (meta-)stable at late times, it may account for the 
dark matter in the universe while allowing for successful 
reheating of the SM particles and setting the necessary 
conditions for BBN. 


III. HYBRID MODEL 


terfall fields. One possibility to satisfy this condition and 
simultaneously generate a constant vacuum energy is to 
introduce an additional “driving” superfield, Z , which is 
not charged under the discrete symmetry and is coupled 
to the waterfall sector, along the lines proposed in [36| |. 
We thus consider a superpotential of the form: 

W = |$(Y 2 - X 2 _) + |$(Y 2 - Y 2 ) + - ^(Y l + Y 2 ) 
+ ^ZiXl + Xl - M 2 ) 

+ ^(X+ + A_)Q 2 + ... , (3.1) 

where M is a constant mass scale and we have included a 
coupling between the waterfall superfields and additional 
chiral superfields Q which give their decay products. The 
dots indicate additional terms that may be added, involv¬ 
ing the inflaton and the superfield Z. In particular, if 
the scalar component of the latter has a sufficiently large 
mass, either from superpotential terms, soft masses from 
SUSY breaking in other sectors or non-minimal terms in 
the Kahler potential, its expectation value will be set to 
zero both during and after inflation. The global mini¬ 
mum of the scalar potential will then lie along the real 
direction (A + ) = (X + ) = x/v2, which preserves the dis¬ 
crete symmetry, and the scalar potential relevant for the 
inflationary and post-inflationary dynamics has the usual 
hybrid form: 


A. Basic properties and dynamics 


V(<t>,x) 



(3.2) 


An alternative framework to the one considered in the 
previous section is the supersymmetric hybrid model 0. 
In this scenario, the inflaton decay products need not 
include or interact with the Standard Model degrees of 
freedom, since the additional waterfall sector can be re¬ 
sponsible for reheating after inflation :32f, 33| . 

For the same reasons exposed in the minimal model, 
the symmetry Ci C 'Ll x Si is imposed on the super- 
field containing the inflaton and all the superfields that 
it couples directly to. Hence, for the superpotential to 
be invariant under the action of this group, the inflaton 
is, as before, coupled to a pair of superfields Y± which 
contain the fermions with masses m /, and to a water¬ 
fall sector with a pair of superfields X±. The group C 2 
simultaneously changes $ —> — $ and interchanges the 
superfields Y + o Y_ and X + o Y_. 

The discrete symmetry forbids the linear term in the 
superpotential that is typically considered in SUSY hy¬ 
brid inflation to generate the constant vacuum energy 
driving accelerated expansion. This may nevertheless 
be generated either by a D-term contribution [34j], or 
through a non-vanishing F-term coming from a SUSY 
breaking sector [35]]. In addition, in order to ensure the 
stability of the inflaton at late times, so that it may ac¬ 
count for dark matter, the C 2 symmetry must be pre¬ 
served in the ground state, implying equal vacuum ex¬ 
pectation values for the scalar components of both wa- 


where (f> = V2{<&) is the real inflaton scalar field. We 
recover the usual SUSY hybrid case for re = g/\J2 and 
for simplicity we will consider this parametric regime, 
although our analysis can be extended to the generic case. 

Inflation takes place for amplitudes of the inflaton field 
larger than a critical value, (f> > (j> c = M/y/ 2, such that 
the waterfall field is held at the origin % = 0. As the in¬ 
flaton rolls towards its minimum at 4> = 0, its amplitude 
falls below the critical value and the waterfall field can 
roll to its true vacuum at x = Af, thus ending inflation. 
After that point, both fields start to oscillate around its 
respective minima, triggering the process of reheating the 
universe into a radiation era. 

In this scenario, the inflaton cannot decay into either 
the bosonic or fermionic components in the waterfall sec¬ 
tor due to kinematical blocking. This is easy to see at the 
global minimum, where to/ = m x = gM, but extends to 
all field values. As in the minimal model, the inflaton can 
decay into the Y± fields and the decay will be incomplete 
for to/ > to//2 = gM/2. For simplicity, we assume that 
the scalar components of the Y± fields acquire large soft 
masses from SUSY breaking and focus on the fermionic 
decay channels, noting that the inclusion of both chan¬ 
nels will not change our conclusions significantly. 

The waterfall fields will decay into the Q sector fields 
and, for similar reasons, we include only the fermionic de¬ 
cay channels in this case as well. We assume that these 
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fields are light, eventually leading to the complete decay 
of the waterfall sector and reheating the universe. Note 
that neither the inflaton nor the waterfall field can be 
completely stable in order to reheat the universe, since 
they carry a comparable amount of the energy density 
after inflation. In particular, if the inflaton were com¬ 
pletely stable and behaved as dark matter at all times, 
the decay of the waterfall field would only convert half 
of the total energy density into radiation. Its incomplete 
decay will then reduce the inflaton abundance and hence 
allow for an efficient reheating once the waterfall field de¬ 
cays. For the inflaton to decay incompletely before the 
waterfall field, we require h > h x , which is the parametric 
regime on which we will focus henceforth. 

The evolution equations driving the post-infiationary 
dynamics of the inflaton and waterfall fields, as well as 
their decay products which we assume to quickly ther- 
malize, are then given by: 

0 + 3170 + g 2 x 2 (f> = -r^, (3.3) 

2 

X + 3 Hx + y (X 2 - M 2 + 20 2 )x = ~T x x, (3.4) 

pji + 4Hpfi = T x x 2 +T c f,<j) 2 , (3-5) 

where the decay width of the inflaton is = F + + T_, 
with r ± given by Eq. m , while the decay width of the 
waterfall field is given by: 


r* = 



(3.6) 


where we neglect the Q fermion masses, which we have 
checked numerically to be a good approximation for cou¬ 
plings h x < 0.1 in the parameter range of interest to our 
discussion. 

Eqs. (13.31) . (13.41) and (13.51) . together with the Friedmann 
equation, 


tt 2 _ P<l> + Px + PR _ + I-X 2 + ^x) + PR 

3Mp ~ 3Mp 

(3.7) 

form a complete set of differential equations that 
can be solved numerically given a set of parameters 
(g, M,m,f,h : h x ) and initial conditions. Fig. [3] shows 
an example of the numerical solution for the parameter 
values M = 10 _2 Mp, g = 10 -5 , h = 1, h x = 10 -5 and 
S = 0.02, the latter being defined in Eq. (12.91) . In this 
example we have considered initial conditions such that 
the fields are close to their values at the end of inflation 
and have small velocities, 0(0) = l.OOO10 c , x(0) = 0, 
0(0) = x(0) = —10 ~ 4 gM 2 , with also pr(0) = 0 since any 
pre-inflationary radiation should be exponentially diluted 
by the accelerated expansion. However, we have checked 
numerically that the results do not show a strong depen¬ 
dence on the choice of initial conditions. 

At the beginning of the evolution, both the inflaton 
and the waterfall field mimic a pair of coupled matter 
fluids which give a roughly equal contribution to the total 
energy density. After the incomplete decay of the inflaton 




FIG. 3. Results of the numerical integration of the inflaton- 
waterfall-radiation dynamical equations for M = 10~ 2 Mp, 
g = ICE 5 , h = 1, h x = 10~ 5 and 5 = 0.02, showing the time 
evolution of (a) the inflaton (solid blue curve), the waterfall 
(dashed green curve) and radiation (dotted red curve) energy 
densities; (b) the inflaton-to-entropy ratio. All quantities are 
given in Planck units such that Mp = 1. 


becomes efficient, its energy density is transferred into 
the radiation bath, while the waterfall field evolves as 
an effective non-interacting matter field in an expanding 
universe. 

Since radiation is diluted more quickly than matter by 
the cosmological expansion, the universe experiences an 
era of matter domination until the waterfall field effec¬ 
tively decays and reheats the universe. We then enter the 
standard radiation-dominated era of Big Bang cosmol¬ 
ogy, with the oscillating inflaton field remnant behaving 
as a cold dark matter component. 

The temperature after the effective decay of the in¬ 
flaton, Tjj, can be determined analytically following the 
same reasoning used in minimal model to compute the 
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reheating temperature, giving: 


T d 



1/4 


( 

\Mp 


3/4 

M P , 


(3.8) 


where g*p> is the effective number of light degrees of free¬ 
dom when the inflaton effectively decays. Unlike in the 
minimal model, the reheating temperature is not con¬ 
trolled by the inflaton but rather by the waterfall decay, 
corresponding to the temperature for which T x = H: 

Tr ~ 0.23g~ft 4 h x , (3.9) 

with g*R being the effective number of light degrees of 
freedom at reheating and where we used = m x since 
the fields are close to the global minimum at this stage. 

The computation of the inflaton dark matter abun¬ 
dance after reheating is more involved than in the min¬ 
imal scenario due to the larger set of parameters in the 
hybrid scenario and the coupling between the inflaton 
and waterfall fields. The effect of the parameters of the 
potential in Eq. (13721) . g and M = m^/g on the final 
abundance is shown in Fig. [2 


the endpoint of the curves at small g values such that 
= yfgM = m^/yfg = 10 16 GeV, of the order of the 
GUT scale, which gives roughly the upper bound on the 
scale of inflation. The vacuum energy scale therefore di¬ 
minishes when moving from left to right along each curve 
of constant in Fig. 2) 

We will restrict the remainder of our analysis to the 
values of g for which the incomplete decay exhibits its 
maximal efficiency, therefore yielding the lowest inflaton 
abundance for a given mass m^. Numerically, we find 
these values to be given by: 

9 ~ °-° 4 ( Sp ) r °' 08 ' (3 - 10) 

With this relation, we obtain for the inflaton-to- 
entropy ratio from the numerical results: 

^ ~ 134 g- 1/A h x h~ 3A2 ^S 2 ^ ( ^V’ 7 ,(3.11) 

s V Mp / 

P(h) = 1 - ^ log 10 h , (3.12) 

j(S) = 1 + O.1151og 10 <5. (3.13) 



FIG. 4. Dependence of the inflaton-to-entropy ratio after re¬ 
heating multiplied by the inflaton mass in Planck units for 
different values of the parameters in the potential, g, M for 
5 = 0.02, h = 1 and h x = 1CT 5 . 

For large values of the coupling g , the inflaton-to- 
entropy ratio grows with g since the system is strongly 
coupled and the waterfall field transfers a significant part 
of its energy to the inflaton, therefore increasing its abun¬ 
dance. For very low values of the coupling, on the other 
hand, the inflaton-to-entropy ratio shows a very mild 
dependence on this coupling. In this case the Hubble 
parameter H ~ m|/g is very large, so that the condi¬ 
tion for the incomplete decay to be efficient, > H, 
cannot be maintained for sufficiently long and the abun¬ 
dance is not so drastically reduced. Numerically, we take 


This expression is equivalent to Eq. (12.171) in the mini¬ 
mal model, with the different powers on the parameters 
reflecting the more complicated dynamics present in this 
scenario. Equating (12.191) and (13.111) . we then obtain 
for the inflaton mass yielding the observed dark matter 
abundance: 


m 0 ~ 442 5 ° 15 °' 6 h -°mh) d -o.<ry(6) GeV 

(3.14) 

For couplings h ~ 0(1) and h x ~ 10 3 — 10 5 , the in¬ 
flaton may account for the dark matter in the universe 
with masses in the GeV-TeV range, similarly to the mini¬ 
mal model, while predicting a reheating temperature well 
above the BBN constraint. The qualitative dependence 
of the inflaton mass on the Yukawa coupling h can be 
understood using the same arguments as in the minimal 
realization for inflaton dark matter described earlier in 
this work. 

The parameter h x determines the duration of the re¬ 
heating process, which is due to the decay of the waterfall 
field. During that process, there is entropy production 
[37j | and as a consequence the abundance n^/s will be 
further reduced, as shown in Fig. [3] (b). The dilution 
factor is given by T = S eq /S r, with S being the entropy 
and the subscript “eg” denoting the time at which the 
X thermalized decay products start dominating the ra¬ 
diation bath, at a temperature T eq > Tr. The dilution 
factor is then given by: 



5 /Tr\ / g*p \ 
3 \Td) \g*R/ 


(3.15) 


with Tp> given in Eq. 13.81 The smaller the coupling 
h x , the longer the reheating process and the smaller Tr, 
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and the more efficient the reduction of the abundance. 
Therefore larger inflaton masses are allowed to match 
the present dark matter abundance. 


B. Condensate evaporation: the WIMPlaton 
scenario 


The effect of additional processes induced by the 
Yukawa coupling h in the evaporation of the conden¬ 
sate of zero-momentum inflaton particles is essentially 
the same that has been described for the minimal model, 
with the rate of evaporation given by Eq. (12.231) . How¬ 
ever, the possibility of the fermions becoming non- 
relativistic during the matter era, where the dominant 
contribution to the energy density comes from the wa¬ 
terfall field, gives a different lower bound on the Yukawa 
couplings such that evaporation is inefficient for T > 
m^, to /: 


h x < 5 x 10- 10 /i-w 1 °sio' 1 


^g*_y 3/32 


(3.16) 


Nevertheless, as we found in the minimal scenario of infla¬ 
ton dark matter, this bound is in tension with the lower 
bound on the couplings consistent with a reheating tem¬ 
perature above 100 MeV. Therefore, in the region where 
the condensate may account for the present dark mat¬ 
ter and the reheating temperature is consistent with the 
BBN constraint, condensate evaporation is most likely 
inevitable. The evaporation process produces a bath 
of fermions and inflaton particles kept in local thermal 
equilibrium by annihilation and elastic scatterings that 
eventually become inefficient, at which point the inflaton 
abundance freezes out as in the standard WIMP scenario. 

If freeze-out occurs after the waterfall field has de¬ 
cayed, i.e. in the radiation-dominated era, we obtain 
the same value for the inflaton mass as in the minimal 
model, given in Eq. (12.281) . Freeze-out takes place in the 
waterfall-dominated matter for values of h x satisfying: 


(g*F \ 3 / 3 , 

pFj" 11 / 8 / 

^ ^00^0 \ 

V 10 ) 

l 25 / ' 

v o.i ) 


1/4 


(3.17) 


which then yields an upper bound on the WIMPlaton 
mass: 


m <£ ^5 756 h 2 



(3.18) 


where the annihilation cross section is given by Eq. 
(12.271) . In Fig. [5] we summarize on the plane (h,h x ) 
the different possibilities for inflaton dark matter in the 
hybrid model. In both regions where the freeze-out of 
the inflaton abundance occurs either in the radiation or 
waterfall era, the inflaton can account for the present 
dark matter abundance for masses in the GeV-TeV range 
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FIG. 5. Parameter space of the hybrid model of inflaton dark 
matter with <5 = 1. In the orange (blue) region the abun¬ 
dance of the inflaton accounts for the present dark matter 
energy density and the freeze-out occurs in the radiation (wa- 
terfall=matter) era. The purple (green) region is excluded 
because the reheating (freeze-out) temperature is below 100 
MeV. Dashed (dashed-dot) lines are curves of constant infla¬ 
ton mass (reheating temperature). The grey area represents 
the transition between the regions where the freeze-out takes 
place in the radiation and waterfall era. 


with the reheating and freeze-out temperatures being 
well above the limit imposed by BBN. This shows that 
the WIMPlaton scenario introduced earlier is not an ex¬ 
clusive feature of the minimal model, with a single dy¬ 
namical field, but also occurs in other models of inflation 
with additional dynamical fields. 

While the inflaton mass values corresponding to the 
observed dark matter abundance are not very different in 
the two realizations that we have analyzed, in the hybrid 
scenario there are novel phenomenological possibilities. 
In particular, the inflaton decay products need not inter¬ 
act with the SM degrees of freedom, since these may be 
excited only after the decay of the waterfall field. Either 
the waterfall sector decays directly into SM particles or 
its decay products interact with some of the latter. We 
note that the waterfall fields may be charged under gauge 
symmetries, in which case the relevant terms in the su¬ 
perpotential are of the form &X±X±, etc, where X± and 
X± transform in conjugate representations of the gauge 
group. This will then open up new avenues for model¬ 
building in inflaton dark matter scenarios besides those 
described in the minimal model, which lie, however, out¬ 
side the scope of this work. 


C. Embedding in an inflationary model 

We have explored in the previous sections the dynam¬ 
ics of the incomplete decay of the inflaton field within 
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a (longer) reheating period controlled by a distinct wa¬ 
terfall sector. The difference with respect to the mini¬ 
mal model is the additional entropy production at the 
end of reheating which dilutes the inflaton-dark matter 
abundance, therefore allowing for larger inflaton masses. 
However, as in the minimal case, evaporation of the con¬ 
densate seems unavoidable, leaving us with similar con¬ 
straints for the inflaton mass as in the previous scenario, 
i.e. masses in the GeV-TeV range. 

As for the minimal model, such a low mass for the in¬ 
flaton field at the minimum constrains the form of the 
scalar potential that simultaneously yields a consistent 
inflationary model and accounts for the observed dark 
matter abundance. In particular, either inflation can oc¬ 
cur at a low energy scale or there must be a large hi¬ 
erarchy between the values of the inflaton mass during 
inflation and at the minimum of the potential. 

In the standard SUSY hybrid models with minimal 
Kahler potential, inflation is driven essentially by the 
constant vacuum energy Vo = g 2 M 4 /8 while the waterfall 
fields are stabilized at the origin. We may then use the 
normalization of the scalar curvature power spectrum: 


pi/2 = JL _L 

K 27 tMp ’ 


(3.19) 


and the observational value P% ~ 5 x 10 5 [23| to obtain 
the relation: 


> 2.5 x 10 15 HGeV , (3.20) 

for 4> > M, where m<£ = gM is the inflaton mass at 
the minimum. Since the scalar spectral index n s — 
1 + 2g ~ 0.9603 ± 0.0073 at 68%CL [23| in these scenar¬ 
ios, we conclude that < TeV cannot yield an obser- 
vationally consistent model. Typically we have from the 
normalization of the spectrum M ~ 10 13 — 10 16 GeV, 
and then the WIMPlaton scenario requires small cou¬ 
plings g ~ 10 -13 — 10“ 10 , which are responsible for the 
very flat potential during inflation and the scale invari¬ 
ant spectrum. A non-minimal Kahler potential can yield 
the observed spectral index for lower values of the cou¬ 
pling, although at the expense of a slight increase in M 
[38l . [39| , which again makes it difficult to achieve the re¬ 
quired WIMPlaton mass values. 

An interesting possibility that can realize inflation at 
low energy scales is to inflate along the waterfall trajec¬ 
tory, which has been shown to yield a sufficiently red- 
tilted spectrum for small couplings g within the context 
of both SUSY and non-SUSY hybrid models [4f||4l|. Al¬ 
though it has been recently realized that isocurvature 
contributions yield a too large enhancement of the am¬ 
plitude of the spectrum unless N e < 60, the WIMPlaton 
scenario typically favors a long reheating period, with 
e.g. N e ~ 30 — 40 being allowed by the data [42| . 

As for the minimal model, a simple way to realize the 
WIMPlaton scenario is to inflate in a chaotic rather than 
in a vacuum-dominated regime [4314451 ] . In particular, a 
quartic self-coupling A 2 0 4 /4 can be easily introduced by a 


superpotential coupling between the inflaton and the the 
driving field, \<& 2 Z/2 as in Eq. (12.3411 . We then have: 

£)’( 

(3.21) 

so that the quartic term will easily dominate over the 
vacuum term for super-planckian values and the typical 
parameters required by the normalization of the spec¬ 
trum and the WIMPlaton scenario. We note that the 
vacuum term will come to dominate the energy at small 
field values, but for (j> > 4> c we have 
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(3.22) 


such that slow-roll inflation never takes place in the small 
field regime. Inflation may then occur entirely in a 
chaotic regime, with a non-minimal coupling to gravity 
or warm inflation yielding observables within the Planck 
window as described for the minimal model. 

We note that, in supergravity models, such chaotic in¬ 
flation scenarios can be obtained by considering a non- 
minimal Kahler potential for the inflaton, while tak¬ 
ing the canonical one for the other superfields in the 
model, in particular the Z field. One possibility is to 
consider a Kahler potential with a shift symmetry [46[ . 
e.g. AT(<f>, Z, ...) = ($ + 4>+) 2 /2 + ZZ^ + ..., with infla¬ 
tion taking place along the imaginary component of the 
scalar inflaton. 


IV. CONCLUSION 

In this work we have shown that the decay of the infla¬ 
ton following the inflationary slow-roll regime can be in¬ 
complete, such that successful reheating is achieved while 
leaving a stable remnant that can account for the ob¬ 
served dark matter in the universe. This is achieved by 
coupling the inflaton field (f> to a pair of fermions ip± 
(and/or similarly scalar fields) while imposing a discrete 
symmetry that simultaneously changes <f> —> — <fi and in¬ 
terchanges the two fermions, corresponding to the dis¬ 
crete subgroup C 2 C Z 2 x S 2 . This symmetry forbids 
all inflaton decay channels except for </> —> ij}±ip± if if is 
preserved in the vacuum state. The inflaton is thus sta¬ 
ble at late times if the bare fermion mass mj > m^j 2, 
where denotes the inflaton mass at the origin. How¬ 
ever, since the physical fermion mass is field-dependent, 
inflaton decay may occur during the initial oscillations 
while the field amplitude is sufficiently large. 

We have then shown that this partial decay is suffi¬ 
cient for a thermalized bath of radiation, which includes 
the fermions and potentially other degrees of freedom, to 
become the dominant component. Since the oscillating 
inflaton condensate does not decay completely, it can ac¬ 
count for the observed dark matter in the universe for 
masses parametrically within the GeV-TeV range in sce¬ 
narios with a single dynamical scalar field. However, by 
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estimating the scattering rate of zero-momentum infla- 
ton particles off thermalized fermions, we concluded that 
the oscillating condensate will most likely evaporate in 
parametric regimes where the reheating temperature is 
above the threshold required for Big Bang Nucleosynthe¬ 
sis. The stable inflaton particles then reach a thermalized 
state which eventually decouples from the cosmological 
radiation bath and freezes out as a standard WIMP. 

In this WIMPlaton scenario, inflaton masses must also 
lie in the GeV-TeV range to account for the observed 
dark matter abundance. While these mass values may a 
priori seem too low to yield the correct amplitude for the 
primordial spectrum of curvature perturbations, we have 
shown that the inflaton mass can be much larger during 
the slow-roll period than at the minimum of the poten¬ 
tial. This is for example the case of large-field chaotic 
models, where the inflaton mass during inflation is set by 
a quartic self-interaction at super-planckian values and a 
bare quadratic term close to the origin. Such models 
are consistent with the latest CMB anisotropy measure¬ 
ments by Planck if e.g. non-minimal couplings to gravity 
or dissipative effects are included. 

In the simplest models with a single dynamical field, 
the inflaton decay products must interact with the Stan¬ 
dard Model (SM) degrees of freedom in order to excite 
them in the thermal bath. We have explored differ¬ 
ent possibilities for such interactions, including the case 
where the inflaton decays into right-handed neutrinos 
or mili-charged particles in a hidden sector. We have 
also explored the alternative possibility of hybrid infla¬ 
tion models, where a dynamical waterfall sector, which 
is also charged under the discrete symmetry, is respon¬ 
sible for reheating the universe. While the incomplete 
decay of the inflaton is still required to allow the rela¬ 
tivistic decay products of the waterfall field to dominate 
the energy balance, in such scenarios the inflaton decay 
products need not couple directly to the SM fields, which 
opens up new phenomenological possibilities. 

Since the entropy produced by the waterfall decay di¬ 
lutes the inflaton condensate’s abundance, it could in 
principle allow for larger masses accounting for the cor¬ 
rect dark matter abundance. However, condensate evap¬ 
oration will also most likely occur in the viable para¬ 
metric regimes, such that a WIMPlaton scenario with 
masses in the GeV-TeV range is again the most proba¬ 
ble outcome in hybrid models. A chaotic inflation regime 
could also provide a consistent inflationary model in this 
case, although there exist viable scenarios for low-scale 


inflation such as along the waterfall direction. 

A possible alternative that may generically allow for 
small inflaton masses is the curvaton scenario [13, HU , 
where the fluctuations of an additional field that is light 
during inflation dominate the total curvature perturba¬ 
tions when decaying. In the standard realization of this 
scenario, the curvaton is a late decaying field, and indeed 
this would bring some additional dilution of the inflaton 
abundance due to entropy production. While this may 
relax the bound on the inflaton mass in the condensate 
scenario, it would have practically no influence when the 
condensate evaporates, and the same limits on the WIM- 
PLaton mass apply. 

Our results show that there exist concrete particle 
physics models where a single dynamical field can drive 
inflation in the early universe and account for cold dark 
matter at late times. This is an appealing feature from 
the theoretical perspective, since it allows one to address 
two of the most important problems in modern cosmol¬ 
ogy within the same simple model. Moreover, we have 
found that consistent cosmological scenarios typically re¬ 
quire the masses of both the inflaton and the particles it 
interacts directly with to lie in the GeV-TeV range, open¬ 
ing up the possibility of testing both inflation and dark 
matter at present particle colliders such as the LHC, as 
well as direct or indirect dark matter searches. Finally, 
this scenario also singles out particular classes of infla¬ 
tionary models that can be tested with CMB experiments 
such as Planck, which are now reaching unprecedented 
levels of precision. We thus hope that our work motivates 
further exploration of the phenomenological and observa¬ 
tional consequences of the proposed generic mechanism 
in its several possible implementations. 
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